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Abstract 


An injectivity test consists of continuously injecting a phase (water or gas), for a period of 
time, into an oil-saturated reservoir. From the analysis of the wellbore pressure behavior, this 
procedure aims to estimate, in addition to the volume of recoverable oil, reservoir parameters, 
such as permeability, skin effect among others. In this context, this work aims to propose an 
approximate analytical solution for the pressure behavior during a water injectivity test on 
a multilayer reservoir considering a scheme of multiple injection flow-rates. The accuracy of 
the proposed solution was evaluated by comparison to a commercial finite difference-based 
flow simulator in different scenarios. Results express a considerable agreement between 
the data provided by the numerical simulator and the proposed model. In addition, from 
the proposed model it was possible to estimate the reservoir equivalent permeability with 
satisfactory results. 


Keywords: Injectivity Test, Multilayer Reservoir, Multiple Flow-Rate, Permeability 
Estimation 





1. Introduction and Background 


Injectivity test is a method based on 
injecting a fluid into a reservoir (usually, 
water) to optimize oil extraction or obtain 
information about the reservoir. From the 
pressure response received in this process, 
relevant information regarding the reservoir 
can be obtained, such as_ equivalent 
permeability, recoverable oil volume, among 
others. 
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Considering a single-phase flow, in 1987, 
Ehlig-Economides et al. [1] propose a 
new method to estimate each reservoir 
layer’s properties and interpret them. Such 
methods are based on the analysis of 
pressure and flow-rate for each layer. Thus, 
it is possible to determine the permeability 
and formation damage of each layer. In 
reservoirs with cross-flow formations, the 
vertical permeability between the layers 
could also be estimated 


In 1989, Raghavan [2] synthesized the 
advances until then regarding multilayer 
reservoirs under single-phase flow. His 
main objective was to extend the 
available pressure analysis methods of 
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commingled fluid production, examining 
the characteristics of standard semilog 
techniques. 

Thompson and Reynolds [3] present in 
1997 a pressure-transient behavior study 
of single and multiphase flow in radially 
heterogeneous reservoirs the light of a 
mobility averaging theory that is based on 
Darcy’s law. Therefore, the model can 
apply to injectivity tests. 

In 2001, Peres and Reynolds [4] propose 
an approximated analytical solution for 
injectivity tests on horizontal wells. The 
results are compared with a numerical 
simulator. Then, it is shown that 
the approximated solution is sufficiently 
accurate for practical applications. Besides 
that, it is shown that one can easily 
construct dimensionless solutions without 
resorting to numerical simulation. 

Approximate analytical solutions for the 
injection pressure response at a horizontal 
water injection well have been presented by 
Peres et al. [5] in 2004. 

Subsequently, Peres et al. [6] used the 
superposition principle to provide a solution 
during the falloff period. Although the 
injection problem is not linear, the authors 
followed Abbaszadeh and Kamal [7] and 
Bratvold and Horne [8] works. They 
demonstrated that the analytical solution 
for the falloff period could be obtained with 
reasonable precision under the assumption 
that the total mobility profile remains 
stationary after the well is shut-in and 
using an expression for the flow-rate profile 
obtained. 

Barreto et al. [9] in 2011 propose 
an analytical solution for vertical water 
injector well completed in a multilayer 
commingled reservoir. This solution is used 
to compute the wellbore pressure, injection 
flow-rates, and the waterfront for each layer. 


Also, parameter estimation can be obtained 
from the proposed solution. 

In 2019, Bela et al. [10] propose an 
analytical solution to the falloff period in 
multilayered reservoirs based on the use 
of a dimensionless flow-rate. The model 
efficiency was verified by comparison to a 
numerical simulator. 

Finally, in 2020 Bonafé et al. [11] 
presented a new formulation to approximate 
the pressure response at the well when 
performing injectivity tests with multiple 
flow-rates in a single layer reservoir. The 
proposed formulation is also used to 
determine the reservoir permeability of any 
specific injection or falloff period. The 
accuracy of the proposed solution was also 
assessed by comparison to a numerical 
simulator. 

This paper proposes an approximate 
solution for the pressure response in 
a multilayer reservoir with a vertical 
well during an_ injectivity test with 
multiple flow-rate schemes. The proposed 
formulation was reached by extending and 
generalizing the existing injection-falloff 
solutions presented by [9], [10], and [11], 
considering the total volume of injected 
water into the reservoir during the test. 
The proposed solution’s accuracy was 
evaluated by comparison to a commercial 
flow simulator, considering a selected set 
of synthetic cases. The model will also 
be used to compute parameters related to 
both water and oil phases for each injection 
period. 


2. Mathematical Model 


Starting from the formulation for an 
injectivity test in a single layer reservoir 
with constant flow proposed by Peres et. 
al [5], Bonafé [11] proposed a formulation 


for an injectivity test in a single layer 
reservoir using a multiple flow-rate scheme. 
Following these works, this paper will 
develop the formulation that describes the 
pressure behavior during water injection in 
a multilayer reservoir considering a multiple 
flow-rate scheme. In this formulation, a 
reservoir with n-layers will be considered. 
Water injection occurs with a multiple flow 
scheme considering m-injection periods. 
The Figure 1 shows an example of this 
reservoir model. 

For the mathematical formulation, the 
following hypotheses are assumed: 


e Reservoir is homogeneous and 
isotropic, with infinite extension; 


e Negligible gravitational and capillary 
forces; 


e Water and oil are considered 
immiscible and _ slightly compressible 
fluids with constant viscosity (1); 


e Isothermal flow; 
e Constant thickness in each layer (h;); 
e The well fully penetrates all layers; 


e The well injects a constant flow rate q; 
at each injection period 1; 


e The rock formation has a small and 
constant compressibility; 


e There is no wellbore storage effect; 


e There is no damaged region close to the 
well (skin); 


e At time ¢t = 0, the reservoir is in 
equilibrium, that is, the pressure is the 
same in all layers. 


Therefore, the flow-rate behavior in the 
well is described according to the following 
scheme: 
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where q;; represents the injection flow-rate 
in layer 7 at the injection period i, while q 
represents the injection flow-rate in the well in 
the injection period i. 


3. Formulation for 
Reservoirs 


Multilayer 


This section synthetize the solution of the 
pressure behavior on a multilayer reservoir 
during water injectivity test found in previous 
works. The fluid mobility (Ay) measures 
the easiness of a fluid to flow through the 
porous medium. So, following the work [3], 
this property is defined as the ratio between 
the relative permeability of the phase and its 
viscosity, that is: 


ky ¢ (Sw) 


A (Sw) = - 


(2) 


where, f = w,o. 
Thus, the total mobility of the oil-water 
system can be defined as: 


At (Sw) = Ao( Sw) a Aw(Sw) (3) 


Using Darcy’s law, oil and water flow-rates 
are given by: 





do, (r,t) = —2hjkjo, (« 2pilr ) (4) 


and 
duy(tot) = —2aghjde, (rE) 6) 


r 


By mass conservation, total flow in layer 7 is 
obtained by: 


qj (r, t) = do; (r, t) + dw; (r, t) (6) 


Thereby, replacing equations (4) and (5) in 
Equation (6): 


Op; (r,t 
OAT) == 2 hgh; (Ao, + Az) (ree?) 
(7) 
During an injection period a_ water 


saturation region is formed around the 
well. Over time, this region propagates, so 
saturation is a function of time and radius. 
Therefore, Equation (7) can be rewritten as 
follows: 

qj(r,t) = —2rhjkj rz, (7, t) (Ee) (8) 


; 
Rearranging the terms of the Equation (8): 


Opj(r,t) L(t) 1 
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From the assumptions the reservoir is 
infinite. Therefore, integrating both sides of 
Equation (9) from the well: 

[o-e) 
1 qj(r,t) dr 
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On the other hand, the water injected 
volume (vj;) at layer j, for tj-1 <t < t;, where 
i =1,...,m, can be defined as: 


wyilt) = f atta (11) 
0 


Thus, considering the Buckley-Leveret 
formulation [12], the waterfront radius (rp,) 


can be defined as: 
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It is important to note that fi,,(Sw) 
represents the water fractional flow derivative 
in layer 7, at each saturation point (S,,). The 
fractional flow curve depends on the mobility 
of the fluid: 


fw; (Sw) = ee (13) 


In each layer, it is easy to notice that in the 
region where r > rp, only oil flows. Thus, 
using the constants of adequacy to the system 
of measures ap, the integral of Equation (9) 
can be reformulated as the following sum: 
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where a,=19.03. 
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Equation (15) shows that the pressure 
variation in the well can be understood as the 
sum of two terms: one due to differences in 
mobility between oil and water (Ap),) and 
another related to the single-phase oil flow 


(Apo; ): 


Apw,(t) = Ap, (t) + Apo, (t) (16) 


where: 
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The approximated solution of Equation (16) 
was presented by [9]. 


4. Proposed Formulation 


To determine the pressure variation along 
the reservoir, an approximation for the 
multiple flow-rate scheme in each layer (q;) will 
be used. Assuming that the flow through each 
layer can be approximated by the flow in the 
flooded region or by the oil flow ahead of the 
flooded region: 


G(r, t) © ay(r,t) © do;(1; t) (19) 


Therefore, the term referring to the mobility 
difference between water and oil in layer 4, 
explained at the previous section by Equation 
(14), can be rewritten as follows: 
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By hypothesis, pressure variation at the well 
in each layer is the same, that is: 


Apus(t) = Apup, (t) =---= Apwfr(t) (21) 


As an algebraic device, the pressure variations 
in the well in each layer can be added. 
This mathematical device, despite having no 
physical sense, helps to obtain a general 
solution for pressure variation at the well. 


From this, using Equations (16) and (21), two 
equalities are obtained: 


p (22) 
= S° (Apy, (t) + Apo, (t)) 
j=l 
So: 
Apus(t) = = | So Apa,(t) + > Apo, (t) 
=1 j=1 
J J (23) 


On the other hand, the single-phase solution 
ensures that Apo, is in fact the same in all 
layers [9]. That is: 


Ap) = Ap (= = Die (24) 
Thus, 
S- Apo, = nApo(t) (25) 
j=l 


Where Ap,(t) is the pressure change 
obtained from the solution for single-phase flow 
problem given by: 
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Finally, replacing Equation (25) in Equation 
(23) it is possible to obtain an approximate 
expression for the pressure behavior in the well: 


a Apy, (t 


Apws(t )+Ap,(t) (27) 


It is important to note that for a single 
layer problem the solution proposed in the 
Equation (27) reduces to the solution proposed 
by Bonafé et al [11]. 


5. Flow-rate scheme approach 


In the formulation proposed in Section 4 
a crucial step is relative to use an approach 
for the flow-rate at each reservoir layer. To 
determine this approach a same reasoning 
proposed by [11] was followed. 

For a single-layer reservoir and considering 
a single-phase flow with constant flow-rate, 
from the combination of the mass conservation 
equation, Darcy’s law and state equations, the 
solution for the pressure is given by: 


AngB 1 rhe 
Ap(r,t) = | i ( - )| 28 
Prt) hkXo 2 AarkXot (28) 





Thus, for a single-phase flow and considering 
a single-layer reservoir with a multiple 
flow-rate scheme, it is possible to apply 
flow-rate superposition. Therefore, the 
pressure solution is now given by: 


Ap(r,t) = 


hag Pos, (rete 
re * ( nat) ) 
29 





On the other hand, by Darcy’s law: 
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So, replacing the Equation (29) in (30): 
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Therefore, in the formulation proposed in 
Section 4, for a multilayer reservoir the 
flow-rate in layer j in a given injection period 
4 can be approximated by profile flow-rate 
obtained for a single phase flow: 


a r2 bce 


BO) => Ne ase ree 
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6. Results and Discussion 


The accuracy of the proposed solution was 
evaluated by comparison with a commercial 
flow simulator based on finite differences. 
Scenarios considering different reservoir 
properties were used in such comparisons. 
To correctly interpret the results of an 
injectivity test, it is necessary to analyze the 
pressure data and the pressure derivative curve 
concerning the natural logarithm of time. For 
numerical calculations, the pressure derivative 
was calculated as proposed by Bourdet [13] 
for each time variation, and for numerical 
integration, the trapezoidal rule was used. 


6.1. Comparison with the Numerical 


Simulator 


For all cases of this section, a total injection 
period of 72 hours (3 days) was considered, 
where each period refers to 24 hours (1 day) 
of injection, according to Table 1. In each case 
group, scenarios with different permeability in 
each layer and differents oil viscosity values 
were considered. The relative permeability 
curves showed in Figure 2 are evaluated for 
each layer. 

Cases A and B were divided into four cases, 
according to Table 2. In these cases, two 
pairs of cases with different reservoir properties 
were considered. The difference between these 
sub-cases is the oil viscosity used, which 
influences the oil mobility value (A,). The 
mobility ratio is a useful variable to interpret 
which fluid moves more easily in the porous 
medium. This parameter has no physical 
meaning and comes from the mobility ratio 
curve M, which expresses the ratio between two 
corresponding points on the mobility curves: 


MS») = Sa (33) 





Thus, the mobility ratio (MZ) will be defined 
as the ratio between the water mobility in the 


residual oil saturation and the oil mobility in 
the initial water saturation: 





MS (34) 


Because of that, in cases where M > 1, 
the flow is favorable to water, while in cases 
where Mf < 1 flow is unfavorable to water. 
The graphs presented in this section are on the 
log-log scale to facilitate the visualization of the 
transient pressure pulse. The presented results 
were divided into two groups of cases according 
to the injection flow-rates schemes. 

Cases A: In this set of cases, the injection 
rate increases gradually at each period. The 
first two cases to be presented, called Cases Al 
and A2, have the same characteristics shown 
in Table 2. For case Al, it is possible to 
observe good agreement between the curves 
in all periods shown in Figure 3. Figure 4 
show the results for case A2, it is possible 
to observe a small discrepancy between the 
derivative curves at a given time. This occurs 
because of the spatial discretization of the flow 
simulator and due to shocks in the solution 
inherent to the problem. Thus, the response 
of the flow simulator is felt later in comparison 
to the analytical response. 

Cases B: In this set of cases, reservoir 
properties are the same as those used in case 
A, but the flow-rate scheme is different. As 
shown in Table 2, injection flow-rates are 
gradually reduced during injection periods. 
Consequently, the pressure curves for the 2nd 
and 3rd periods are decreasing. In general, 
although the injection flow rates in Cases B are 
different from those in Cases A, the behaviors 
and the causes of the curves’ errors are similar. 
It is also possible to observe a good similarity 
between the curves in Figures 5 and 6 for Cases 
B1 and B2, respectively. 

In general, the graphs showed a good 
correspondence between the curves of the 
derivatives of the model built in this work and 
the curves generated by the flow simulator. 


Thus, it is possible to conclude that the model 
presented for the problem is useful. The 
formulation presented can be used for modeling 
a large number of flow-rates or reservoir layers. 


6.2. Estimated Equivalent Permeability 


The proposed solution can also be 
applied to estimate the reservoir equivalent 
permeability.To evaluate the accuracy of 
the estimated permeability, two long-term 
injectivity tests with 516 hours-flow were 
considered. The cases in this section are 
called C and D and use the multiple flow-rate 
scheme in Table 3. Besides that, the reservoir 
properties used as input parameters can be 
seen in the Table 4. 

The reservoir equivalent permeability is 
calculated using the permeabilities of each 
layer inserted in the simulation. Thus, the 
reservoir equivalent permeability is given by 
the following expression: 


kyhy +...+knhn 
fae 7 





(35) 
On the other hand, for the single-phase flow 


case with constant flow-rate, it is known that 
the pressure change at the wellbore is given by: 


ApgB AarkXot 
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It is easy to notice that Equation (36) 
represents a linear graph against the natural 
logarithm of time whose inclination is given by: 


apqB 


Mo = 1.151 (37) 


O 

The line slope allows estimating the reservoir 
equivalent permeability for single-phase flow. 

Adapting this interpretation to a multiple 
flow-rate scheme, the change of injection 
flow-rate (Agq;) will be used. This is defined 
as the absolute value of the injection flow 
difference of two subsequent periods: 


Age = |9:= G1) (38) 


Graphs of pressure data normalized by 
the change of the injection flow-rate help to 
identify the similarity between the pressure 
derivative levels. Such levels, can be associated 
through an interpretation method with the oil 
or water properties in Figures 7, 8, 9, 10. 

Similarly to the interpretation for single 
phase flow, it is possible to estimate the 
reservoir equivalent permeability. Analyzing 
the pressure logarithmic derivative curve at 
each injection period, the properties of each 
fluid are used. Figures 7, 8, 9, 10 show the 
normalized pressure and derivative pressure 
curves for all periods. 

It is important to note that pressure 
behavior in the well is similar to the 
single-phase oil flow during the first injection 
period since the region close to the well is 
saturated with oil. However, in subsequent 
periods, this region is now saturated with 
water. Thus, the well bottom-hole pressure 
behavior becomes similar to the single-phase 
water flow. 

Therefore, in general, for two-phase flow, 
two slopes corresponding to the oil and water 
mobilities must be considered, depending on 
the injection period. Because of that, for cases 
C1, C2, D1, and D2, Figures 11, 12, 18 and 
14, respectively, show the slopes correspondent 
to water and oil properties from the derivative 
curve. With these values, it is possible to 
calculate the estimated reservoir permeability 
using the relative permeability (water and oil) 
endpoints. 


u AqB 
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Thus, the results obtained in cases 


C1, C2, D1 and D2 for the calculated 
equivalent permeability and the percentage 
error compared to the reservoir equivalent 
permeability can be seen in Tables 5, 6, 7 and 
8, respectively. 


The results obtained for estimated 
reservoir permeability compared to reference 
permeability, in general, are very close for all 
periods. These results prove the effectiveness 
of the proposed model to estimate the reservoir 
equivalent permeability. 


7. Conclusions 


Based on the formulation for the falloff 
period by [10] in multilayer reservoirs, and 
the formulation for injectivity tests with 
multiple flow-rates by [11], it was possible to 
propose a solution for the pressure behavior 
on multilayer reservoirs considering multiple 
flow-rate schemes in this work. The developed 
analytical model can also be used as a basis 
for future work. As an example, the presented 
solution can be extended to a solution in the 
Laplace Domain for multilayer reservoirs and 
multiple flow-rate schemes. Moreover, the 
model presented can also be used as a basis 
for similar formulations considering a damage 
region (skin) or temperature effects. The 
proposed model was applied in a set of cases 
with different reservoir properties, such as a 
different number of layers, permeabilities, and 
injection flow-rates. The comparison between 
the approximate analytical solution and the 
numerical simulator indicated good agreement 
for all tested cases. Furthermore, it was also 
possible to estimate the reservoir equivalent 
permeability during each injection period. The 
low error obtained for the estimated equivalent 
permeability proves that the proposed model 
effectively estimates this reservoir parameter. 


Nomenclature 


ce = Total compressibility 

fw = Fractionary flow 

f. = Fractionary flow derivative 

hj(m) = Thickness of layer j 

ht(m) = Reservoir total thickness 

keq(mD) = Reservoir equivalent permeability 
kj(mD) = Permeability in layer j 


M = Endpoint mobility ratio 

Pw (kgf/cm?) = Wellbottom hole pressure 
ging(m?/day) = Injection flow-rate 
qj((m?/day)) = Injection flow-rate in each 
layer j 

rpj(m) = Waterfront radius in layer j 
Tw(m) = Wellbore radius 

Sw = Water saturation 

t(h) = Time 

Q, = Pressure unit conversion constant 

az = Time unit conversion constant 

= Euler constant 


No = Endpoint oil mobility 

Aw = Endpoint water mobility 
At = Total mobility 

o = Porosity 

[lo(cP) = Oil Viscosity 
Lw(cP) = Water Viscosity 
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Figure 1: Reservoir model with n layers and two-phase flow 


0.6 


0.5 


Relative Permeability 
—) —) 
Nh w 


0.1 





0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Water Saturation (S,) 


Figure 2: Relative permeability curves 
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Figure 4: Pressure and pressure derivative curves for case A2 
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Figure 5: Pressure and pressure derivative curves for case B1 
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Figure 6: Pressure and pressure derivative curves for case B2 
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Figure 7: Normalized pressure and pressure derivative curves for case Cl 
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Table 1: Flow-rate schemes for cases A and B 
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Figure 8: Normalized pressure and pressure derivative curves for case C2 
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Table 2: Rock and fluid properties for cases A and B 
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Figure 9: Normalized pressure and pressure derivative curves for case D1 
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Table 3: Flow-rate schemes for cases C and D 
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Figure 10: Normalized pressure and pressure derivative curves for case D2 
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Table 4: Rock and fluid properties for cases C and D 
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Figure 11: Normalized pressure derivative curves for case Cl 
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Figure 12: Normalized pressure derivative curve for case C2 
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Table 5: Estimated equivalent permeability results for case C1 
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Table 6: Estimated equivalent permeability results for case C2 
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Reg ; Keg Error Bist Error 

Case (mD) Period (mD) (%) (mD) (%) 
1 749.90 -0.01 741.002 -1.20 

2 798.26 6.43 717.31 -4.36 

3 902.41 20.42 718.13 -4.25 

PE} 8 | 4 | 796.08 | 614 | 720.23 | -3.97 
5 Tiel 3.04 726.55 3.13 

6 815.12 8.68 735.57 1.92 























Table 7: Estimated equivalent permeability results for case D1 











Keq Reos Error ae Error 
Case (mD) Period (mD) (%) (mD) (%) 

1 749.30 -0.09 745.74 -0.57 

2 752.04 0.27 744.34 -0.75 

3 773.01 3.07 748.60 -0.19 

nel os 4 705.54 -5.93 747.14 -0.48 
5 672.76 -10.40 745.73 -0.57 

6 756.30 0.84 742.78 -0.96 
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Table 8: Estimated equivalent permeability results for case D2 








